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We consider the problem of the crossover from BCS superconductivity to Bose-Einstein condensation 
in three dimensions for a system of fermions with an attractive interaction, for which we adopt the 
simplifying assumption of a suitably regularized point-contact interaction. We examine in a critical 
way the fermionic (self-consistent) T-matrix approximation which has been widely utilized in the 
literature to describe this crossover above the superconducting critical temperature, and show that 
it fails to yield the correct behaviour of the system in the strong-coupling limit, where composite 
bosons form as tightly bound fermion pairs. We then set up the correct approximation for a "dilute" 
system of composite bosons and show that an entire new class of diagrams has to be considered in 
the place of the fermionic T-matrix approximation for the self-energy. This new class of diagrams 
correctly describes both the weak- and strong-coupling limits, and consequently results into an 
improved interpolation scheme for the intermediate (crossover) region. In this context, we provide 
also a systematic mapping between the corresponding diagrammatic theories for the composite 
bosons and the constituent fermions. As a preliminary result to demonstrate the numerical effect 
of our new class of diagrams on physical quantities, we calculate the value of the scattering length 
for composite bosons in the strong-coupling limit and show that it is considerably modified with 
respect to the result obtained within the self-consistent fermionic T-matrix approximation. 



PACS numbers: 74.25.-q, 74.20.-z, 05.30.Jp 

I. INTRODUCTION 

The problem of the crossover from BCS superconduc- 
tivity to Bose-Einstein (BEl condensation has attracted 
considerable interest lately,B and especially after the re- 
cent ARPES experiments in cuprate superconductors 
which have shown the existence of a (pseudo) gap at tem- 
peratures-above the superconducting transition temper- 
ature T c .EfH This observation has prompted the proposal 
by many authors of the possible presence of quasi-bound 
fermionic pairs above T c and up to a second temperature 
scale T* . More generally, one may think of ascribing the 
anomalous experimental properties of the cuprate super- 
conductors in the normal state (i.e., for temperatures T 
such that T c ~ T ~ T*) as well as the possible break- 
down of the FermiJiquid theory to the presence of these 
quasi-bound pairs .□ 

From the theoretical point of view, this crossover prob- 
lem (as well as related crossover problems, like the one 
associated with the Mott transitions) poses a compelling 
challenge, because approximations that are valid on the 
one side of the crossover are not necessarily valid on the 
opposite side. In addition, the crossover region is char- 
acterized (or even denned) by the absence of a "small" 
parameter, which would allow one to control the approx- 
imations. For systems that are sufficiently "dilute" (such 
that, for instance, for given strength of the inter-particle 
interaction, the density p can be taken to be arbitrarily 
small), one could directly exploit the well-known results 
obtained for "dilute" fermionic systems with pp = p, 



on the one hand,Eil or for "dilute" bosonic systems with 
Pb = p/2, on the other handj3 to get a correct description 
of the limits on the two sides of the crossover. 

For a "dilute" fermionic system the two-body equation 
plays an essential role.Q In particular, in three dimen- 
sions the low-energy two-body scattering process can be 
parametrized in terms of the scattering length ap, which 
is negative for weak coupling and positive for strong 
coupling (i.e., in the presence of a bound state for the 
attractive interaction), and diverges when the coupling 
strength suffices for the bound state to appear. For 
strong coupling, of gives the size of the bound state. The 
many-body approach to a "dilute" fermionic system then 
introduces the small dimensionless parameter kpap to 
select the relevant diagrammatic contributions, whereby 
the Fermi wave vector identifies the average interparticle 
distance kp 1 and ap emerges directly from the many- 
particle T-matrix (which is defined by repeated scatter- 
ing in the particle-particle channel at finite fermionic den- 
sity). In the original treatment with a purely repulsive 
interaction, no bound state occurred and clf remained 
finite, in such a way that the parameter kpap could un- 
ambiguously be taken to be much less than unity.13 In the 
case of interest to us of an attractive interaction, which 
develops a bound state in three dimensions, however, the 
parameter A^ofI may readily exceed unity when the 
strength of the attractive interaction increases toward 
its critical value for the appearance of a bound state. 
In this case, the many-body T-matrix can no longer be 
identified with ap. Furthermore, past the critical inter- 
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action strength where the two-body problem develops a 
bound state (with ap turning again finite and eventually 
vanishing in the strong-coupling limit), the many-body 
T-matrix acquires a singularity (pole) for all interaction 
strengths and no identification of the T-matrix with ap is 
clearly possible any longer. Therefore, although the pa- 
rameter fcj? aj? may still serve to readily locate which side 
of the crossover one is examining (or even qualitatively, 
how close to the crossover region one is) by relying on the 
results of the two-body problem, the smallness of this pa- 
rameter might not by itself be sufficient to guarantee that 
an approximation, selected for the weak-coupling limit, 
is still valid in the strong-coupling limit. 

In particular, the T-matrix approximation for the 
fcrmionic self-energy (which has invariably been regarded 
in the literature as representing the dilute-approach ap- 
proximation to the BCS-BE crossover problem)0lI3 is not 
expected to remain valid as soon as the two-body bound 
state develops, since for this approximation the power 
counting in the small parameter kpap relies just on the 
identification of the many-body T-matrix with ap. This 
statement is consistent with the physical picture that, 
as soon as the two-body bound state develops, it is the 
residual interaction between the composite bosons to de- 
termine the "diluteness" condition of the system and not 
the original attraction between the constituent fermions, 
which produces the bound state to begin with and rel- 
atively to which the system is in the strong-coupling 
limit £.3 It is then clear that a correct description of the 
strong-coupling limit can be obtained only by selecting 
the relevant approximations directly for a dilute system 
of composite bosons, rather than relying on the fermionic 
T-matrix approximation, which is valid by construction 
for a dilute system of fermions. 

From previous work on the functiopalpititegral ap- 
proach to the crossover from BCS to BE,E-jH3 one knows 
that approximations (like the BCS mean field at zero 
temperature), which give a satisfactory account of the 
weak-coupling limit, become inadequate in the strong- 
coupling limit, and that only by including fluctuation cor- 
rections at least at the one-loop level a sensible descrip- 
tion of the effective bosonic system in the strong-coupling 
limit results. This remark has actually-suggested to deal 
with the crossover problem in reverse Jl3 that is, by first 
envisaging approximations which give a satisfactory de- 
scription of the strong-coupling (bosonic) limit and by 
extrapolating them toward the weak-coupling (fermionicl 
limit, where they are expected to work properly as well.tS 

One then anticipates the set of many-body (self- 
energy) diagrams, which describe the bosonic limit, to 
be much richer than the corresponding set of diagrams, 
which describe the fcrmionic limit. In fact, from the pre- 
vious discussion we do not expect the fermionic T-matrix 
approximation to be appropriate for describing the "di- 
lute" bosonic limit properly. We will indeed show below 
that the fermionic T-matrix approximation corresponds 
in the strong-coupling limit to the standard bosonic 
Hartree-Fock approximation, and that, for this reason, 



it misses all but one of the infinite set of self-eneijgy dia- 
grams associated with a "dilute" bosonic systcmEj 

The dynamics of a "dilute" system of true (point- 
like) bosons can be accounted for by the correspond- 
ing T-matrix approximation (aside from a narrow re- 
gion about T c ).LlEH Using a representation with a sym- 
metrized two-body interaction, the diagrams correspond- 
ing to the bosonic T-matrix approximation contain a se- 
quence of pairs of propagators running in the same di- 
rection, with only one additional propagator running in 
the opposite direction (see Fig. ||b of subsection III-A be- 
low). The first term of this sequence (i.e., the one with no 
pairs but only with the single propagator running back- 
ward) corresponds to the standard bosonic Hartree-Fock 
approximation.!^ Since all other terms (with an arbitrary 
number of pairs running in the same direction) are of 

the same order in the "gas parameter" p^as (where 
as is the (positive) scattering length for the low-energy 
two-boson scattering problem) , keeping the Hartree-Fock 
term only is evidently not justified for a "dilute" bosonic 
system. 

On physical grounds, for a "dilute" system of com- 
posite bosons (obtained by pairing fermions via the at- 
tractive interaction in the strong-coupling limit) one ex- 
pects the same picture to emerge (as it will also be con- 
firmed by the results of the present paper). For this sys- 
tem, the bosonic single-particle propagator corresponds 
to a two-fermion Green's function in the particle-particle 
channel (with opposite spins for a spinless boson), the 
analogue of the bosonic self-energy corrections being as- 
sociated in the fermionic particle-particle channel with 
insertions of the many-body effective interaction. The 
complete correspondence between the two (bosonic and 
fermionic) diagrammatic structures requires one to estab- 
lish a well-defined mapping between the building blocks 
of the two structures (namely, propagators and inter- 
action vertices). One of the purposes of this paper is 
to discuss this mapping between the two diagrammatic 
structures, thus complementing the mapping established 
by means of functional integrals in Ref. 16. It is worth 
mentioning in this respect that the correspondence we 
shall find between the symmetry factors of the bosonic 
diagrammatic structure, on the one hand, and the num- 
ber of independent diagrams in the associated fcrmionic 
structure, on the other hand, provides a compelling check 
on our mapping procedure. 

The reason for choosing the particle-particle fermionic 
channel in order to single out the (approximate) dy- 
namics of the system in the strong-coupling limit stems 
from the obvious fact that this channel only has physi- 
cal relevance in this limit. The single-particle fermionic 
propagator (or, equivalently, the fermionic irreducible 
self-energy) is then obtained as a derived quantity, by 
suitably closing the two-fermion diagrams with a sin- 
gle fermionic . line. In this way, the Baym-Kadanoff 
procedureE3t3 is effectively inverted, since we construct 
the one-particle (fermionic) self-energy diagrams in or- 
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der to recover the chosen set of two-particle diagrams 
via functional differentiation, rather than following the 
standard procedure of deriving the two-particle diagrams 
from a chosen set of one- particle diagrams O 

As already mentioned, til the interaction between the 
composite bosons effectively vanishes in the extreme 
strong-coupling limit. In this extreme limit, neither 
the Hartree-Fock bosonic correction (which results from 
the standard fermionic T-matrix) nor the set of "low- 
density" bosonic self-energy diagrams (which are in- 
cluded in the present approach) need then to be taken 
into account, as they both correspond to interaction cor- 
rections. The many-body corrections introduced in this 
paper in the bosonic limit are thus expected to be im- 
portant whenever the interaction between the composite 
bosons has direct influence on physical properties, such 
as the critical tempetature T c and the bosonic binding 
energy (pseudogap).H3 In addition, including the inter- 
action between the composite bosons should make our 
"low-density" approximation to be valid somewhat closer 
to the crossover region, over and above the boundary set 
by retaining the Hartree-Fock term only. In this respect, 
our finding that the composite-boson scattering lenght 
as is smaller than the value = 2of (obtained in Ref. 
9 within the fermionic T-matrix approximation) effec- 
tively restricts the range of the crossover region on the 
BE side, and makes the eventual extrapolation of our ap- 
proximation through the crossover region more reliable. 

In this paper, we present the formal theory for the 
choice of the fermionic self-energy diagrams starting from 
the strong-coupling side, for temperatures above T c (in 
the sense that no anomalous single-particle fermionic 
propagator will be considered) and for three spatial di- 
mensions. The corresponding theory below T c remains 
to be developed. An extensive numerical study based on 
this approximation is under way and will be discussed 
separately. The only numerical calculation presented in 
this paper concerns the value of the composite-boson 
scattering length as in the strong-coupling limit. 

The plan of the paper is as follows. In Section II we dis- 
cuss some introductory material, which is necessary for 
setting up our "low-density" approximation for compos- 
ite bosons. Specifically: We introduce a suitable regular- 
ization for the fermionic point-contact interaction, which 
allows us to select readily the relevant classes of fermionic 
diagrams; We summarize the mapping onto a bosonic 
system in the strong-coupling limit, obtained by the pro- 
cedure of Ref. 16; We discuss the standard fermionic T- 
matrix approximation and show how the Hartree-Fock 
approximation for composite bosons results from it in 
the strong-coupling limit. In Section III, after recall- 
ing the standard theory of the "low-density" Bose gas 
for true (point-like) bosons, we introduce the theory of 
the "low-density" Bose gas for composite bosons and 
obtain the fermionic self-energy from the bosonic self- 
energy in the "low-density" approximation. In Section IV 
we present the numerical calculation for the composite- 
boson scattering length as in the strong-coupling limit 



and discuss the physical implications of our result. Sec- 
tion V gives our conclusions. The Appendices discuss 
more technical material. Specifically, in Appendix A we 
prove that, for our choice of a "contact" fermionic in- 
teraction, the composite-boson propagator coincides ex- 
actly with the fermionic generalized particle-particle lad- 
der. In Appendix B we discuss in detail the mapping 
between the bosonic and fermionic diagrammatic struc- 
tures and demonstrate how the symmetry factors of the 
bosonic theory translates onto the number of indepen- 
dent diagrams in the fermionic theory, by working out 
two examples in detail. 

II. BUILDING BLOCKS OF THE 
DIAGRAMMATIC STRUCTURE FOR 
COMPOSITE BOSONS 

In this Section we discuss the diagrammatic structure 
that generically describes the composite bosons in terms 
of the constituent fermions, as a preliminary step for set- 
ting up the "low-density" approximation for composite 
bosons in the next Section. Our construction rests on a 
judicious choice of the fermionic interaction, which (al- 
beit without loss of generality) greatly reduces the num- 
ber and considerably simplifies the expressions of the 
Feynman diagrams to be taken into account. We will 
show how this construction fits with the general map- 
ping established in Ref. 16 via functional-integral meth- 
ods. We shall, in turn, utilize the results of Ref. 16 to 
single out the lowest-order (four-point) bosonic vertex as 
the relevant one in the strong-coupling limit. We shall 
further show that the standard fernairmic T-matrix ap- 
proximation for a "dilute" Fermi gasOu reproduces in the 
strong-coupling limit the Hartree-Fock approximation for 
composite bosonsjla which can be defined in analogy, with 
the standard result for true (point-like) bosons .EjO In 
this way, we shall make manifest the necessity of super- 
seding the standard fermionic T-matrix approximation 
when dealing with the crossover problem and of replac- 
ing it by a more complete approximation, which enables 
us instead to recover the "low-density" approximation for 
composite bosons in the strong-coupling limit. 

A. Regularization of the fermionic interaction 

We begin by considering the following simple model 
Hamiltonian for interacting fermions (we set Planck h 
and Boltzmann ks constants equal to unity throughout): 

+ ^E/ drdr'^(r)^(r')^eff(r-r / )^(r')^(r) (2-1) 

a, a' 
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where tp a (r) is the fermionic field operator with spin pro- 
jection a = (|, I), m the fermionic (effective) mass, \i the 
fermionic chemical potential, and V c s(r — r') an effective 
potential that provides the attraction between fermions. 

To simplify the ensuing many-body diagrammatic 
structure considerably (and yet preserving the physical 
effects we are after) , we adopt for V^ff the simple form of 



a contact potential 



Kff(r-r') = v J(r-r') 



(2.2) 



where v is a negative constant. With this choice, the in- 
teraction affec ts o nly fermions with opposite spins in the 
Hamiltonian (2.1) owing to P auli principle. A suitable 
regularization of the potential (|2.2| ) is, however, required 
to get accurate control of the many-body diagrammatic 
structure. In particular, the equation (in the center-of- 
mass frame) 
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(2.3) 



for the fermionic scattering length of associated with the 
potential (2.2) is ill-defined, since the integral over the 
three-dimensional wave vecto r k is ultraviolet divergent. 
The delta- function potential (2.2) is then effectively reg- 
ularized, by intr odu cing an ultraviolet cutoff k Q in the 
integral of Eq. (2.2) and letting v Q — > as k Q — > oo, 
in order to keep Of fixed at a chosen finite value. The 
required relation between v a and k Q is obtained directly 
from Eq. (tO). O ne nnds: 



where f3 is the inverse temperature, w„ = (2n + l)7r/3 _1 
(n integer) a fermionic Matsubara frequency, Q° a (bare) 
single-particle fermionic Green's function, fp(E) = 
{e 0E + ly 1 the Fermi function, f(k) = k 2 /(2m) - fj,, 
and R pp is a remainder that converges to a finite value in 
the limit k — > oo. Note that I pp = — I + 0(k~ 1 ) owing 



to Eq. (2.4) 



k,+k 



k 2 -k 



(a) 



k' 



kj+k ' 
— s — 1 — >- 



f k,+k-k' 
i 



k,-k' 



k 2 -k 



k,+k 



k,+k 



(b) 
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when fc |aF| 3> 1. Note that a? enters the expression 
(2.4) only through the subleading term proportional to 
k~ l . Keeping this subleading term turns out to be es- 
sential for our purposes. [Note that, if v D were positive 
instead of negative, the scattering length of would van- 
ish in the limit k Q — * oo even keeping v finite.] 

With the regularization ( p.4| ) for the potential, the clas- 
sification of the many-body diagrams gets considerably 
simplified, since only specific sub-structures of these di- 
agrams survive when the limit k Q — > oo is eventually 
taken. Let us consider, in particular, the characteristkj 
rungs occurring in the theory of a "dilute" Fermi gas,B 
which are depicted in Fig. |l|a and Fig. [j]b, respectively. 
For the particle-particle rung of Fig. |l|a we obtain: 



Ipp(ki,w ni ; k2, lu. 



-y 



elk 



v(k) 



x g°(k 1 +K,uj ni 



P^J (2tt)3 
■f k, uj ni + u) n ) G°(k 2 - k, uj n2 - u n ) 

k °dk [M^ki + k)) + f F (gk 2 - k)) - i] 

(2tt) 3 i(w Ba + u ni ) - C(ki + k) - e(k 2 - k) 
mk 



2tt 2 



w rai ;k 2 ,o; n2 ) 



(2.5) 



(c) 



(d) 



FIG. 1. (a) Particle-particle and (b) particle-hole rung; 
(c) vertex correction and (d) two-particle effective interac- 
tion. Full and broken lines represent the "bare" single-particle 
fermionic Green's function and the fermionic interaction, re- 
spectively. Four-momenta are indicated and spin labels are 
represented by up and down arrows. 

By the same token, for the particle-hole rung of Fig. lb 
we obtain: 

1 f elk 

7 ph .(ki,cj ni ;k2,cj„ 2 ) = J j^v^) (2.6) 

x g°(ki + k, oj ni + u) n ) G°(k 2 + k, uj n2 + u n ) 

• £ °dk [/Hflki+k)) - f F (gk 2 + k))) 
(2nY i(uj n2 - u ni ) + C(ki + k) - ^(k 2 + k) 

which vanishes in the limit k —f oo owing to the vanish- 
ing of the prefactor v . In this case, the integral over k on 
the right-hand side of Eq. (2.6) remains finite in the limit, 
because fp (£(k)) is proportional to exp {— (3k 2 /(2m)} 
for large |k| (irrespective of the value of the chemical 
potential). This situation has to be contrasted with the 
original argument by Galitskii, whereby for a finite-range 
potential the particle-hole rung I p h does not vanish, but 
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it is reduced with respect to its particle-particle coun- 
terpart J pp by a factor kpap [hp being the Fermi wave 
vector) £3 In the weak-coupling limit, in fact, the integral 
on the right-hand side of Eq. (2.6) is proportional to kp 
from dimensional considerations. 

In an analogous way, one can show that in the particle- 
particle channel the contributions of the vertex correc- 
tions (like the one depicted in Fig. |l|c) and of the two- 
particle effective interactions other than the rung (like 
the one depicted in Fig. |l]d) vanish for our choice of 
the potential, since they both contain a factor I p h of the 
particle-hole type. [Note that the diagram of Fig. con- 
tains also a forbidden interaction between parallel spins.] 

Despite these drastic simplifications, care should be ex- 
erted in not dismissing too promptly contributions which, 
albeit being individually vanishing as k Q — > oo, belong to 
an infinite sequence that could instead converge to a finite 
value. Two examples of alternative behaviors are given 
in Fig. ^a and Fig. 0b. In particular, for the particle- 
particle ladder of FigTga we obtain from Eq. ( |2.5| ) (apart 
from an overall minus sign originating from the fermionic 
diagrammatic rules) : 



and k a , respectively. One then obtains the following fi- 
nite result for the particle-particle ladder (2.7) 



1 



1 + 0{k F a F ) + 0{k a a F 



(2.10) 



which permits one to classify the diagrammatic structure 
in powers of kpap (and/or of k a ap) in the weak-coupling 
limit J 
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(2.7) 



since both numerator and denominator in Eq. ( |2.7| ) van- 
ish like k^ 1 when k Q — > oo. [Note that we have intro- 
duced the four- vector notation k = (k, uj n ). Note also 
that we have attached the | spin to the upper line and 
the I spin to the lower line of the ladder, respectively. 
This convention will be consistently mantained in the 
following.] For the particle-hole bubble series of Fig. gb 
we obtain instead a vanishing result, since 



1 + lph(q\0) 







(2.8) 



3jS k — 

The particle-particle ladder (2/7), upon surviving the 
limit k a — > oo, represents a building block of the many- 
body diagrammatic structure. It is thus relevant to ex- 
amine in detail its analytic behavior in the weak- and 
strong-coupling limits, in the order. Let 



E,, 



(ki+k 2 ) s 

4m 



+ 2/i (2.9) 



(b) 

FIG. 2. (a) Particle-particle ladder and (b) series of parti- 
cle-hole bubbles. Note that for a point-contact potential the 
particle-particle ladder depends only on the sum of the incom- 
ing (outgoing) four-momenta. Conventions are as in Fig. |lj 

In the strong-coupling limit (whereby Pfj, — > — oo),@ 
on t he o ther hand, only the —1 term within brackets in 
Eq. (2.5) contributes to R pp , yielding 



,3/2 



4-7T 



/ £«sgn(ImE Q ) (2.11) 



with the standard convention for the cut of the square 
root along the negative real axis. Recalling further that 
2|/i| = (map) -1 in the (extreme) strong-coupling limit ,ta 
we obtain for the particle-particle ladder (|]7]) the follow- 
ing polar structure 
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Stt 
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+ R pp (h +k 2 ) m 2 a F i(uj ni +uj n2 ) 



(ki 



4m 



be the (complex) energy entering the denominator in 
Eq. ( |2.5| ) and k a — y / 2mE a the associated (complex) 
wave vector. In the weak-coupling limitj23 bot h th e Jf 
terms and the —1 term within brackets in Eq. (2.5) con- 
tribute to R ppi yielding contributions proportional to kp quency f2„ = L0 ni 



(2.12) 



which (apart from the residue being different from unity) 
resembles a "free" boson propagator with Matsubara fre- 
j n2 — 2vnf3~ 1 (y integer), wave 
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vector q = ki + k2, mass 2m, and vanishing chemical 
potential. [More generally, the bosonic chemical poten- 
tial [Ib = 2u + e should be added to the denominator 



of Eq. ( 2.12 ), where e G is the bound-state energy of the 
two-fermion problem - cf. subsection IIIB. In the extreme 
strong-coupling limit \xb vanishes.] The complete anal- 
ogy between the fermionic generalized particle-par ticle 
ladder (of which the "bare" particle-particle ladder (2.7) 
is only a piece) and the "full" boson propagator is dis- 
cussed in Appendix A. 

It is evident from the above considerations that the 
classification of diagrams developed for a "dilute" system 
in weak coupling (which relics on the finite value ( 2.10| ) 
of the particle-particle ladder - cf. Ref. 6) can no longer 
be utilized in the strong-coupling limit, since in this case 
the particle-particle ladder does not reduc e to a constant 
but develops a polar structure [cf. Eq. ( 2.12J )]. In the 
strong- coupling limit, a different classification scheme is 
therefore required to organize the diagrammatic structure 
for a "dilute" system. 





FIG. 3. Four- and six-point vertices for composite bosons. 
Incoming and outgoing bosonic four-momenta are indicated, 
as well as four-momenta on each fermionic line. Spin labels 
are understood to alternate on successive fermionic lines. 

It is also evident from the above considerations that, 
with our choice of the fermionic interaction, the skeleton 
structure of the diagrammatic theory can be constructed 
only with the particle-particle ladder (2.7) plus an infinite 
number of interaction vertices, like the ones depicted in 
Fig. |^ (besides one spare single-particle fermionic Green's 



function that enters the fermionic self-energy diagrams, 
which in turn contains in principle all self-energy inser- 
tions originating from self-consistency). 

The ladder and the vertices of Fig. g contain, by con- 
struction, only "bare" single-particle fermionic Green's 
functions. These vertices serve to connect the ladders 
among themselves, thus generating complex diagram- 
matic structures. In analogy with the so-called Hikami 
vertices occurring in the weak-localization problem,c3 we 
refer to these vertices as the four, six, -point vertices, 
in the order. 

As far as the dressing of the particle-particle ladder 
is concerned, the above statements can be justified in a 
rigorous fashion by resorting, e.g., to functional-integral 
methods. To this end, in the following subsection we 
recall some relevant results obtained in Ref. 16, adapting 
them to the purposes of the present paper. 



B. Mapping onto a bosonic system via functional 
integrals 

We briefly review the procedure of Ref. 16 to extract 
the set of effective bosonic interactions mentioned above 
from the original fermionic action. We first note that 
comparison with the results of Ref. 16 requires us to in- 
terpret the Fouri er t ransform of the effective attractive 
potential in Eq. (2.1) as a "separable" potential, of the 
form 



Vw(k)w(k') 



(2.13) 



where V is the quantization volume, V is a negative con- 
stant, and 



tu(k) 



1 if |k| < k D 
otherwise 



Note that the limit k a — » oo recovers our "contact" po- 
tential, provided we identify V = v /V, with v given by 
©■ 

In Ref. 16 the original fermionic partition function was 
written in the form of a functional integral 

Z = J VcDc exp{-5} (2.14) 

with action 

r0 



S 



(It I ^c CT (k,T)JU CT (k,T)+#(T) I . (2.ir» 

k,cr 



where r is the imaginary time a nd (c , c) are Grassmann 
variables. Owing to the form ( 2.13| ) of the interaction 
potential, the quartic part of the Hamiltonian was further 
expressed in terms of the operator 
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B(q,r) = ^ W (k) Ci (-k + q/2,r)c T (k + q/2,r) 



(2.16) 

With the use of the following Hubbard-Stratonovich 
transformation (in the particle-particle channel) 

exp{^%r)6(q,r)}=-l /"d6*(q,r)d6(q,r) (2.17) 



x exp <j - |6(q, r) | 2 + 6(q, r)B(q, r) + 6* (q, r)S(q, r) 



/?]>> 2 (k-q/2) 

k 

V 



[/ F (£(k))+/ F (£(k-q))-l] 



- ^(k) - e(k 



l + 



q) 

(2.25) 



where 6(q, r + (3) = 6(q, t) for any q and r, the Grass- 
mann variables were eventually integrated out, to obtain 
(apart from irrelevant numerical constants) 



Z = J Vb*Vb exp{-5 cf f} 
with the "effective" bosonic action 

S * = -^5>(9)| 2 - trlnM - 



(2.18) 



(2.19) 



In Eq. (2.19), we have set q = (q, f2„), the matrix M is 
given by 



e(k)S k ,k' 



-e(-fc)4,fc' 

(2.20) 



with e(fc) = iui n — £(k), and the trace is over the four- 
momentum and the matrix indices. The effective action 
S e R contains all powers in the bosonic-like variables b. 
Expanding the logarithm in Eq. ( 2.19| ) in fact yields 



oo 

S cS = -trlnMs - —^Mqtf + ]T -trA 2 



1=1 



(2.21) 



where 



Mm/.-/.') = ( e( J } _ e (_ k) '2.22 i 



and 



X(k,k') = w 




b-(fc-fc') 



(2.23) 



Consider the quadratic terms first. The ensuing 
quadratic action reads: 



5 $ = ^Ei & (?)i 2 ^) 



(2.24) 



where 



with I pp given by Eq. (2.5). The "bare" bosonic propa- 
gator is thu s expressed in terms of the particle-particle 
ladder (2.7), by writing 



< b*(q)b(q) >„ ( 2) 







V 1 



Ipp(q) 



T (q) (2.26) 



where T (q) is the particle-particle ladder depicted in 
Fig. ||a (including now the appropriate overall sign that 
was missing in Eq. (2.7)). This identification is shown 
graphically in Fig. ||a for a definite choice of the spin la- 
bels (this convention will be mantained in the rest of the 
paper). Note that the factor /3/V of Eq. (2.26) enters also 
the definition of the "full" bosonic propagator discussed 
in Appendix A. 

It is important to emphasize that the identification 
( |2.26D (as well as the other results of this subsection 
and of Appendix A) holds irrespective of the value of the 
fermionic scattering length a F . Nonetheless, referring to 
a system of interacting composite bosons acquires phys- 
ical meaning in the strong-coupling limit only. 



For the quartic terms we obtain instead from Eq. (2.21) 



S$ - ,trA 4 - — £ 



4 

x u 2 (qi 



q4)b*(q 1 )b*(q 2 )b(q 3 )b(q i ) 



(2.27) 



where the (four-point) effective two-boson interaction is 
given by [cf. Fig. § 



V\ 2 



Mil • • ■ 94) = <5 gi +92,<?3+94 ( ~p J 

x y 1 

4^ e(-k)e(k + q 2 )e(-k + q x - q 4 )e(k + q 4 ) 



(2.28) 



with the limit k — > 00 already taken in the last ex- 
pression. It is worth noting the following features of the 
abov e exp ressions: (i) The factor 1/4 (instead of 1/2) in 
Eq. ( 2.27 ) applies to a symmetrized bosonic two-body 
interaction, which mak es . th e counting of the b osonic 
Feynman diagrams easierEJEII; (ii) The interaction (2.2J) 
is indeed symmetric under the interchange q\ <-> q 2 or 
93 94 j as it can be explicitly verified; (Hi) The in- 
teraction (2.28) depends on wave vectors as well as on 
Matsubara frequencies, revealing in this way the com- 
posite natur e of t he bos ons; (iv) The factors (3 and V en- 
tering Eqs. Q2.27D and fl2.28j ) are needed to comply with 
the standard bosonic and fermionic diagrammatic rules 
and with the definition of the "full" bosonic propagator 
given in Appendix A; (v) The energy denominators in 
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Eq. (2.28) correspond to single-particle (bare) fermionic 
Green's functions, si nce Q °(k) — e(fc) -1 ; (vi) The factor 
of 2 in the definition (2.28) corresponds to the two differ- 
ent sequences of spin labels that can be attached to the 
four fermionic Green's functions, as shown graphically in 
Fig. [l|b (where the identification with the effective two- 
boson interaction is also indicated) . Keeping track of the 
spin labels will, in fact, prove important in the following 
to establish the desired mapping between the bosonic and 
fermionic diagrammatic structures. 

k+q f k'+q 



shown in Fig. ||. 



<=> 



1 
v 



(a) 



H 4 4, 

X 




(b) 

FIG. 4. Graphical correspondence (a) between the "bare" 
propagator for composite bosons (represented by a thick line) 
and the fermionic particle-particle ladder of Fig. ^a, and (b) 
between the effective two-boson interaction and the four-point 
vertex of Fig. ^, where now the spin labels have been explicitly 
indicated in the internal lines. Note that the fermionic lines 
composing the four-point vertex have been rearranged with 
respect to Fig. in order to resemble the bosonic vertex 
more closely (accordingly, the fermionic lines never intersect 
each other in the four-p oint vertex ). A ppropriate powers of 
P/V as required by Eqs. ( |2.26| ) and (2.28) have been indicated 
explicitly. 

When considering the sum S~g + S^g of the 
quadratic and quartic actions, the bosonic propagator 
< b*(q)b(q) >„ (2) „(4) can be expressed in terms of the 



"bare" bosonic propagator (2.26) and of the effective two- 
boson interaction ( |2.28| ) via Wick's theorem. The topol- 
ogy of the resulting diagrammatic structure, as well as 
the symmetry factor of each diagram, a re ., id entical to 
those obtained for true (point-like) bosons HEJ The asso- 
ciated fermionic diagrammatic structure can then be con- 
structed whenever needed by the correspondence rules 




1 
v 




F«o 


k+q' 
s k+q 


r(q) 


I I 
+ 

T t 


r(q) 


' k+q 

-k k+q' -k 


T(q) 




FIG. 5. Graphical correspondence between the compos- 
ite-boson propagator and the two-fermion Green's function in 
the particle-particle channel, to first order in the four-point 
interaction vertex. This vertex can be identified from Fig. |j| 
by setting q2 = qs, = q and q\ = q± = q therein. 

As an example of this correspondence, we show in 
Fig. ||the bosonic propagator to first order in the interac- 
tion, together with the associated fermionic diagrams for 
the two-fermion Green's function in the particle-particle 
channel. Since the symmetry factor of this bosonic dia- 
gram is unity, only two fermionic diagrams are associated 
with it. [The value of the bosonic symmetry factor and 
the number of independent fermionic diagrams are inti- 
mately connected, as it will be discussed in Appendix B.] 
Note that a factor l/(/3V) is attached to the sum over q' 
in both bosonic and fermionic diagrams. Note also that 
the minus sign, which is associated with one power of the 
interaction in the bosonic diagram, is associated instead 
with the presence of a closed loop in the corresponding 
fermionic diagrams (the minus signs associated with the 
fermionic interaction being already taken into account in 
the definition of the particle-particle ladder). Note fi- 
nally that the bosonic self-energy insertion of Fig. |^ has 
the same topological structure of the bosonic Hartree- 
Fock self-energy diagramJ23'Ell 

A typical value of the two-boson effective interaction is 
obtained by considering the strong-coupling limit f3fi — > 
-co and setting q\ = • • • = = in Eq. ( 2. 28] ). One 
gets: 
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, 2 ,ih = 2 



m 2 ap 
8tt 



u 2 (0) 



where! 







«2(0) = 



4?r(2a F ) 
2m 



(2.29) 



(2.30) 



The factor rr^ap /(87r) in Eq. (2.29) reflects the difference 
between the true bosonic propagator and the particle- 
part icleladder in the strong-coupling limit [cf. also 
Eq. (2.12)]. Owing to this difference, u 2 (0) given by 
Eq. Q2.30| ) (and not u 2 (0) given by Eq. Q2.29| )) has to be 
identified with the boson-boson interaction at zero four- 
momenta. We shall returtuto the difference between u 2 
and u 2 in subsection IIIB.Ea 

Recalling further that the scattering length a^ orn 
within the Born approximation, obtained for a pair of 
true bosons (each of mass 2m) mutually interacting via 
a two-body potential with Fourier transform ^2(0) at 
zero wave vector, is given by a^ orn = 2mzi2(0)/(47r), 
Eq. (2.30) yields the following relation between the 



bosonic and fermionic scattering lengths: 

2a F ■ 



,Born 



(2.31) 



The result (2.31) was also obtained in Ref. 9 within the 
fermionic self-consistent T-matrix approximation (which 
corresponds to the bosonic Hartree-Fock approximation 
of Fig. ^ in the strong-coupling limit - see the next sub- 
section), where it was regarded to be the value of the 
scattering length as for a "dilute" system of com posite 
bosons. We will show in Section IV that the result (2.31) 
actually differs from as, when all diagrams associated 
with a "dilute" system of composite bosons are taken 
into account. 



Beside the four-point vertex (2.28 ), the composite na- 
ture of the bosons produces (an infinite set of) additional 
vertices. Let us consider, in particular, the six-poi nt ver - 
tex defined by keeping the terms with I = 3 in Eq. (2.21). 
We obtain: 



<? (6) - 1 trX e - 1 



E 

qi—qe 



xu 3 (q 1 ...q 6 )b*(q 1 )b*(q 2 )b*(q 3 )b(q 4 )b(q 5 )b(q 6 ) (2.32) 

where the (six-point) effective three-boson interaction is 
given by [cf. Fig. § 

W3(9l • • ■ 96) = <V +92+93,94+95+96 ( 7J ) 7^T7 X/ 



(-1) 



e(—k)e(k + q 2 )e(-k - q 2 + 95) 
1 

e(k + q 2 + q 3 - qz)e(-k + qi - qi)e(k + q±) 



(2.33) 



In the strong-coupling limit, whereby is the relevant 
energy scale in the problem, from dimensional considera- 
tions we get u 2 (P/V) 2 ~ |/ir 3/2 and u 3 {(3/V) 3 ~ M" 7 / 2 



(in three dimensions). F or th is reason, the contribution 



of the tree-boson vertex ( 2.33 ) is suppressed with respect 
to the contribution of the two-boson vertex (2.28). To be 
more precise, one should compare the values of similar di- 
agrams constructed with the four- and six-point vertices, 
respectively (like, for instance, the ones depicted in Fig. |5| 
and Fig. ^J). The value of the diagram of Fig. ||is smaller 
than the value of the diagram of Fig. |^ by the quantity 



H- 7/2 4 



- 3 / 2 e|. /2 |mI 1/2 



{k F a F f 



(2.34) 



Here, the factors containing the Fermi energy ep [= 
k F /(2m) — (37r 2 p) 2 / 3 /(2m)] originate from the bosonic 
cycles (cf. Section III), while the fact ors \n\ and |/x 



1/2 



originate from the residue in Eq. ( 2.12 ). The diagram of 
Fig. |^ can thus be neglected in comparison to the dia- 
gram of Fig. ^[ since kpa F <C 1 in the strong-coupling 
limit. 




1 
v 













t 








T 


r(q) 


k+q 


k+q' 
-k-q+q" 


k+q 




i 




V T ) 




i 



FIG. 6. Graphical correspondence between the compos- 
ite-boson propagator and the two-fermion Green's function 
in the particle-particle channel, to first order in the six-point 
interaction vertex. This vertex can be identified from Fig. ^ 
by setting q 2 = qe — q, qi — qi — q', and q 3 = qs — q" 
therein. 

The above argument can be made more general, by 
showing that all interaction vertices can be neglected 
in comparison-with the four-point vertex in the strong- 
coupling limit!EH In this limit, one can thus construct all 
diagrams representing the two-particle Green's function 
in the particle-particle channel in terms of the "bare" 



9 



ladder and of the four- point interaction vertex only. 
This is precisely what one would expect on physical 
grounds, since the interactions involving more than two 
bodies become progressively less effective as the com- 
posite bosons overlap less when approaching the strong- 
coupling limit. However, it is also clear from the above 
considerations that in the extreme strong-coupling limit 
(whereby olf — > 0) the four-point vertex can be neglected, 
too, and the composite bosons become effectively freeH3 
In the next subsection we will show how the (self- 
consistent) fermionic T-matrix approximation can be ex- 
amined in terms of the four, six, • • •, -vertex functions in 
the strong-coupling limit. 



k' 




I I I I ~ 

I I I I I 

I I I I I 

I I I I I 




+ 



(a) 



1 1 1 1 1 
1 1 1 1 1 



C. Fermionic T-matrix approximation in the 
strong-coupling limit 

The T-matrix approximation for a "dilute" Fermi gas 
represents one of the few cases in the many-body the- 
ory where the choice of the self-energy diagpans can 
be controlled by an externaLsmall parameter!^! In the 
original version by Galitskii,tj the T-matrix approxima- 
tion was conceived for a repulsive fermionic interaction 
of finite range (thus excluding bound states) and with 
the scattering length ap always positive (albeit small). 
The fermionic self-energy diagram associated with this 
approximation is depicted at the left in Fig. 0a, and 
is obtained by closing the particle-particle ladder with 
a single-particle Green's function in the only possible 
way. [The diagram at the right ip Fig. 0a was included 
by Galitskii original treatment ,Ej but vanishes for our 
choice of the attractive potential since it contains for- 
bidden interactions between parallel spins. By the same 
token, no spin summation needs to be considered for 
the fermionic loop at the left in Fig. 0a.] In Fig. 0a 
all single-particle lines are regarded to be self- consistent, 
and thus contain self-energy insertiqas of the same kind 
of the ones depicted in the fi.gure.t2l In this way, the 
self-consistent fermioniG,|X-matrix is "conserving" in the 
Baym-Kadanoff sense.E3E3 Recalling that (with our reg- 
ularization of the potential) the particle-particle ladder 
depends only on the sum of the incoming (outgoing) four- 
momenta, the self-energy of Fig. 0a reads 

£ F (fc) = -JL£)r.(fc + fc') W (2.35) 

where T s is obtained from r o by replacing everywhere 
the "bare" Go with the self-consistent Q. 



(b) 



(d) 



FIG. 7. (a) Self-energy diagrams corresponding to the 
self-consistent fermionic T-matrix approximation (full lines 
here represent self-consistent fermionic single-particle Green's 
functions and spin labels have been suppressed); Self-energy 
corrections entering the particle-particle ladder, obtained by 
contracting the (b) four-point vertex and (c-e) six-point 
vertex (full lines now represent the "bare" fermionic sin- 
gle-particle Green's functions). 

In the approximation of the lowest order, self-energy 
insertions are neglected and the "full" single-particle lines 
of Fig. 0a are replaced by "bar e" ones. The particle- 
particle ladder T s of Eq. (2.35) thus reduces to r o of 
Fig. 0a, which (apart from a sign) is approximately given 
by Eq. ( 2. 10] ) in the weak-coupling limit (irrespective of 
the sign of ap). We thus obtain to this order: 



4 0) w 



iiraF p 
m 2 



ep ( — kpap + 0((kpap) 

37T 



(2.36) 



where the density factor p [= fcj,/(37r 2 )] is supplied by 
the closed fermionic loop. Cons idering self-energy inser- 
tions on top of the result ( 2.36 ) modifies Y^p jtp only to 
0((k F ap) 2 ). 

For the attractive fermionic interaction of interest, the 
scattering length ap changes from being negative when 
the two-body problem fails to support a bound state, to 
being positive when the bound state is eventually sup- 
ported by increasing the interaction strength. As al- 
ready discussed in subsection IIA, the sign change of 
ap has dramatic effects on the "bare" particle-particle 
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ladder r o , that changes from the behavior (2.10) in the 
weak-coupling limit to developing a pol ar structure in 
the strong-coupling limit (cf. Eq. (2.12)). In the lat- 
ter limit, the self-energy acquires the density factor 
from the particle-particle ladder and not from the closed 
fcrmionic loop. This can be explicitly seen by transform- 
ing the Matsubara frequency sum in the expression for 
Y>p into a contour integration in the complex frequency 
plane as follows: 



4 o) (*0 = 



o0Z 



1 



(2.37) 



where C is a contour encircling the imaginary axis in the 
positive sense. In the strong-coupling limit, Q°(k',z — 
iu> n ) has a pole for Re(z) > so th at th e contribu- 
tion of this pole to the integral in Eq. (2.37) is strongly 
reduced by the Bose factor exp{— The dominant 
contribution to the integral in Eq. (2.37) thus comes from 
the pole of r o (k + k', z), yielding 



4 o) (*0 



8tt 



dk' 



1 



'»-"/' J (2^) 3 iWn +£(k')-^^ 
1 



e /3(k+k') 2 /(4m) _ I 

8tt 1 p 

m 2 a,F iijJ n + £(k) 2 



(2.38) 



It is interesting to note that the expression ( [2.38 ) for 
the self-energy leads to a double- fraction structure in the 
single-particle fcrmionic Green's function, since 



G(k,u n ) = 



iuj n - £(k) - 



(2.39) 



K(k) 



with B = 4irp/ (m 2 ap). It is then clear that the expres- 
sion ( 2.39| ) has the form of the BCS Green's function, 
namely, 



g(k,u n ) = 



u 2 (k) 



v 2 (k) 



lU>r> 



E(k) 



E(k) ' 



where 



(2.40) 



(2.41) 



weight v 2 (k)). This secondary peak is relevant to pho- 
toemission experiments and can be interpreted as a 
"shadow-band" structure, associated with a real gap in 
the single-particle spectrum; this gap, in turn, might 
be regarded as anticipating the pseudo-gap actually 
observed in ARPES experimentS|-above the supercon- 
ducting transition temperature.ma It is also for this 
promising feature that several authors have considered 
the fermionic T-matrix approximation as representing 
the "dilute-approach— approximation" for the BCS-BE 
crossover problem.ETtJ 

Self-energy insertions associated with the self- 
consistency of the fermionic single-part icle G reen's func- 
tion modify both r o and Q° in Eq. (2.37). As far as 
the self-energy insertions on r o are concerned, they can 
be interpreted in terms of the four, six, • • •, -point in- 
teraction vertices discussed in the previous subsections. 
Typical examples are shown in Figs, ^b-e, where the 
"bare" single-particle lines associated with the vertices 
are marked by arrows. Note, in particular, that the dia- 
gram of Fig. Mb corresponds to a contraction of the four- 
point vertex, while the diagrams of Figs, ^c-e correspond 
to all possible contractions of the six-point vertex. Addi- 
tional diagrams not shown in the figure would then con- 
tain interaction vertices of higher order. Note also that 
the diagrams of Figs. Rb-c have been already considered 
in Fig. @ and in Fig. pT respectively. From the results of 
the previous subsection we then conclude that only the 
diagram of Fig. [?]b needs to be retained in the strong- 
coupling limit, the other diagrams being suppressed with 
respect to it at least by a factor (fcpcii?) 3 (cf. Eq. (2.34)). 
The diagram of Fig. [7]b corresponds to the Hartree-Fo 
approximation for the self-energy of composite bos ons. 
Note finally that self-energy insertions on Q° in Eq. (2.37) 
could be interpreted similarly, in an open-ended way. 

The above argument leads us to the conclusion that 
the fermionic T-matrix approximation reproduces the 
Hartree-Fock approximation to the self-energy of com- 
posite bosons in the strong-coupling limit . This self- 
energy modifies the pole of T Q in Eq. ( 2.37| ) and affects 
eventually the (pseudo) gap structure in the fermionic 
self-energy. There exist, however, additional contribu- 
tions to the self-energy of composite bosons which are 
of the same order (in the small parameter kpap) of the 
Hartree-Fock approximation just discussed. These con- 
tributions are not included in the fermionic T-matrix ap- 
proximation and must be considered separately, as shown 
in the next Section. 



has a gap given approximately by + B/(2\(i\) and 
where 



v 2 (k) = l-u 2 (k) 



(k F a F ) 



3tt [l + (\k\a F ) 2 } 2 



(2.42) 



Besides the main peak at E(k) (with weight u 2 (k)), 
the s pectr al function associated with the Green's func- 
tion ( 2.4C ) has also a secondary peak at —E(k) (with 



III. T-MATRIX APPROXIMATION FOR 
COMPOSITE BOSONS 

In this Section we shall set up an approximation for the 
fermionic self-energy, which describes the "low-density" 
regime both in the weak-coupling (fermionic) and in the 
strong-coupling (bosonic) limits on equal footing. To 
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this end, we rely on the results of the previous Section 
and consider first the construction of the two-fermion 
Green's function in the bosonic limit, in terms of the 
"bare" particle-particle ladder and of the effective two- 
boson interaction. Before examining this procedure in 
detail, however, it is instructive to briefly recall some 
well-known results concerning the self-.grie.rgy for a "di- 
lute" system of true (point-like) bosonsJj'EJtil and in par- 
ticular the method for selecting the relevant self-energy 
diagrams for this system. In this way, we will learn how 
to generalize this method to the case of a "dilute" system 
of composite bosons, for which the effective two-boson in- 
teraction depends on frequency as well as on wave vector. 

We will eventually extrapolate the approximation thus 
obtained from the bosonic to the (fermionic) weak- 
coupling limit through the crossover region. In the latter 
region (where the parameter kpap exceeds unity), how- 
ever, the "low-density" approximation has to be regarded 
merely as a scheme which interpolates between the two 
controlled limits. 



A. Low-density approximation for true (point-like) 
bosons 

The standard argument to select the diagrams giv- 
ing the leading contribution to the self-enepsj-Jbr a 
"low-density" Bose gas proceeds as follows.Ql230 Let 
^(qii Q2, q3, q.4) be the (symmetrized) bosonic interac- 
tion potential, assumed to be vanishing for |q^| ~ r~ x 
(i = l,---,4), where r is the range of the potential. 
We shall also consider temperatures not too higher than 
the BE critical temperature, so that we shall assume 
T ~ p 2 / 3 . Under these conditions, it turns out th&t,every 
cycle in a diagram contributes a factor T 3 / 2 ~ p£3 

To show this, we consider the generic self-energy dia- 
gram of Fig. ^a, containing one closed path constructed 
with "bare" bosonic propagators. We perform explicitly 
the sum over the four-momentum q running along the 
cycle. Let £s(q) = q 2 /(2ms) — jUg, where tub is the 
bosonic mass and [1b the bosonic chemical potential, re- 
spectively. To the diagram of Fig. ||a we then associate 
the expression 

£(ckxt) = E FiP, Q) E frfa Q) -o —f i \ ( 3J ) 

P,Q q ^ B[(l > 

1 1 

V . . . 

i{Q u £s(q + qi) i(fi„ + CV) - £s(q + q r ) 

where t/(q; Q) stands for the product of all interaction 
potentials attached to the cycle (with all wave vectors 
other than q indicated collectively by Q) and F(P, Q) 
for the expression of the rest of the diagram (appropriate 
factors (/3V) _1 are understood to be associated with the 
four-momentum sums). 




(c) 

FIG. 8. (a) Generic self-energy diagram for a system of 
true bosons, containing one cycle (full lines and full cir- 
cles here stand for the "bare" bosonic propagators and the 
(symmetrized) interaction potential, respectively, while dot- 
ted lines represent connections to the rest of the diagram); 

(b) T- matrix approximation to the self-energy of true bosons; 

(c) Elementary rung for the bosonic T-matrix. 

The sum over the Mats ubara frequency r2„ can be read- 
ily performed in Eq. ( |3.l| ), since the interaction potential 
for point-like bosons is frequency independent. For this 
sum we obtain the result: 

1 1 

" e^(i) - 1 fs(q) + in vi - £s(q + qi) 
1 



£s(q) + i$l Vr - £s(q + q r ) 

1 1 
" e Kb (q+qi) _ 1 £ B ( q + qi ) _ - £ B (q) 



£s(q + qi) + *(fiiv - ^vi) -£s(q + qr) 

When integrating further over the wave vector q in 
Eq. (a nd provided /is is negative), each Bose factor 
in Eq~H introduces a cutoff proportional to T 1 / 2 ~ 
p 1 / 3 about q + q^ ; « (with q^ = 0, qi, • • • , q r ). This 
cutoff is much smaller than the range r^ 1 of the bosonic 
interaction potential, owing to the "low-density" condi- 
tion p x ^r <C 1. We can then approximate q ~ every- 
wherein the integrand associated with the first term of 
Eq. ( |3.2| ), q+qi ~ in t he integrand associated with the 
second term of Eq. (|3.2[), and so on, except for the Bose 
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factors that retain the wave vector dependence indicated 
in Eq. (3.2). We thus obtain the following app roximate 



result for the sum over q in the expression (3.1 
1 U(0;Q) f dq 



2m B 



1 



1 



_3L J (2tt) 3 e«*(q) - 1 

2m B 



-if2„ 



2 HI £-5 



(q r -qi) 2 



1 



s /3?s(q+qi) - l 



(3.3) 



(2tt) 3 

where each integral equals the bosonic density pB = p/2. 
For generic values of the wave vectors qx , • • • , q r and of 
the fre quen cies £l Vl , • • • , £l Vr entering Eq. (3.1), each term 
in Eq. ( 3~3| ) is thus propor tion al to the density. The pres- 
ence of one cycle in Eq. (3.1) makes, in turn, the whole 
expression proportional to p. 

It then follows that, for a "low-density" Bose sys- 
tem, the leading self-energy diagrams contain the min- 
imum number of cycles, like the diagrams shown in 
Fig. ||b, which (pfleetijtutes the so-called bosonic T-matrix 



approximation.ErE3'EJ It is interesting to contrast the be- 
haviour of the cycle (whereby the Bose factor introduces 
a cutoff proportional to T 1 / 2 ~ p 1 ^ 3 ) with the behaviour 
of an elementary rung (drawn in Fig. |^c) in the diagram- 
matic representation of the bosonic T-matrix, with one 
potential vertex and two bosonic propagators running in 
the same direction. In this case, the cutoff originates 
from the potential and not from the Bose factor. Let us, 
in fact, assume, for simplicity, that the bosonic interac- 
tion potential is (approximately) constant (~ u(0)) up 
to the cutoff r~ x . It can then be readily shown that, at 
low enough temperature, the dominant contribution to 
the diagram of Fig . ^|c is given by u(0)r~ 1 , in analogy 
with the result ( |2.5| ) for fermions. No density factor thus 
arises in this case. 

To generalize the above results to a "dilute" system of 
composite bosons, it is essential to take into account the 
frequ ency dependence of the effective two-boson interac- 
tion ( 2.28 ), which ma kes the calculation of the analogue 
of the expression (3.1) somewhat more involved. We will, 
however, show in the next subsection that this quantity 
remains proportional to p. 



B. Low-density approximation for composite bosons 

The generic self-energy diagram containing one cycle 
and associated with the propagator < b*(q)b(q) >s cff 
for composite bosons is shown in Fig. ^a. Here, thick 
lines and squares repr esent the "bare" propagator < 
b*(q)b(q) > <2> of Eq. ( |2.26| ) for composite bosons and 



the effective two-boson interaction given by Eq. ([2.28 ), 
respectively, according to the conventions of Fig. ET[Re- 
call from subsection IIB that only the two-boson interac- 
tion needs to be taken into account in the strong-coupling 
limit. 1 




q+q 



(a) 



+ 



(b) 




(c) 

FIG. 9. (a) Generic self-energy diagram for a system of 
composite bosons, containing one cycle (thick lines and full 
squares stand for the "bare" propagators for composite bosons 
and the effective two-boson interaction, respectively, accord- 
ing to the conventions of Fig. ^, while dotted lines represent 
connections to the rest of the diagram); (b) T-matrix approxi- 
mation to the self-energy of composite bosons; (c) Elementary 
rung for the T-matrix of composite bosons. 



The expression 
propagator < b*{q)b{q) >, 



of the "bare" 

(2) is obtained by inserting in 

that 



Eqs. (|2.26| ) and fl2.7p the expression fl2.1l|) for R pp , 
holds specifically in the strong-coupling limit (whereby 
(3p —>■ -co). Writing y2^sgn(ImE Q ) = i\/— E a and re- 



calling the definition (2.9) for E a , we obtain 





< b*(q)b(q) > (2 ) = -- 

°off V 



-1 



Ana i 



47T 



-in„ + u- - 2/i 



D(q) 



qi 

4m 



(J<B 



(3.4) 



where we have used the definition ps — 2/i + e Q for tba 
bosonic chemical potential and introduced the notationtll 



D(q) 




-*n„ + - p B 



V m ap 



(3.5) 



Replacing further if2„ by the complex frequency z, we 
note that, near the position z = q 2 /(4m) - ^ of the 
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pole in Eq. D(q) of Eq. (3J) reduces to ( -0/V 

times) the constant value 8n / (m ap) entering Eq. ( 2.1 2| ) 
and thus plays the role of a residue. Away from this 
pole, D(q) has a cut along the positive real axis for Re(z) 
> -2(i. 



To retain the structure of Eq. (3.1) as close as possi- 
ble (with tub = 2m in the "bare" bosonic propagators), 
it is convenient to transfer the factors D(q) from the 
"bare" propagators < b*(q)b(q) >„(2> to the effective two- 

boson interactions, which these propagators are joined 
to. Specifically, every propagator < b*{q)b{q) >„<2> as- 

sociatcd with an "internal" line of a diagram tranfers a 
factor ^/D(q) to each of the two interaction vertices it is 
joined to, while each propagator associated with one of 
the two "external" lines transfers one factor \J D(q) to 
the single interaction it is joined to and assigns the re- 
maining factor y/D(q) as a proportionality factor to the 
"full" propagator, which the diagram is meant to rep- 
resent. In this way, the "full" propagator acquires the 
same overall factor D(q) of the "bare" propagat or ( 3.4) 
and the effective two-boson interaction of Eq. ( [2.28 ) is 
multiplied by y/D(q 1 )y/D{q2)y/D(q 3 )y/D{qi). We are 
then led to reseating the effective two-boson interaction 
as follows 

2u 2 (gi ■ ■ ■ <?4) (3.6) 
= U 2 (ll ■ ■ ■ qi)y/D{qi)y/D(q2)y/D(q 3 )y/D( qi ) 



with M 2 (gi . . . (74) given by Eq. ( 2.28| ). [Note that the 
above definition accounts, in particular, for the difference 
between G 2 (0) and "2(0) in Eq. ( |2.29[ ).] 

To the generic self-energy diagram of Fig. ||a, contain- 
ing one cycle of "bare" propagators for com posi te bosons, 
we thus associate an expression similar to (3.1), namely, 



£(q ext ) =£V(P,Q)^E/( 9 ;Q) 
1 



1 



6?(q) 
1 



(3.7) 



+ - £s(q + qi) + fii, r ) - £s(q + q r ) 

where U (q; Q) now stands for the product of all rescaled 
effective two-boson interactions 2u 2 attached to the cycle. 
Although the sum over the Matsubara frequency Q, u in 
Eq. ([3.7]) cannot be performed explicitly as we have done 
in Eq. (JO), it can be conve nientl y transformed into a 
contour integration as in Eq. ( 2.37 ). We thus obtain for 
this frequency sum: 



1 



dz 

c 2ni z - £s(q) z - 
1 



1 



i£i Vl - £s(q + qi) 



X 



z + in Vr -6s(q+qr) e^-i 



(3.8) 



To proceed further, we need to spell out the analytic 
properties of U(q,z;Q). From Fig. ||a one realizes that 



each factor w 2 entering the expressi on for U in Eq. (3.7) 
has the following structure [cf. Eq. ( 2.28Q] : 



U2(q + qj,q a ,qb,q + qj-i)=S qi+ q a , qj _ 1+qb O^j 



E 



■f e(—k)e(k + q a )e(-k + q 3 - g j _ 1 )e(fc + q + qj-i) 



(3-9) 



where q a and q^ belong to the set Q o f Eq. (3.7). In other 
words, the cycle variable q in Eq. (3.9) appears only in the 
last factor e(k + q + qj- 
complex variable z of Eq. ( 3 



1 . Viewed as a function of the 
each factor u 2 has thus 
a cut along the positive real axis for Re(z) > — fi. This 
property, in turn, remains tru e fo r the rescaled effective 
two-boson potential it 2 of Eq. ( |3.6[ ) , since we have already 
determined that the factor D{q) has itself a cut along the 
positive real axis for Re(z) > — 2(i. 

It is thus clear that the contribution of th e singulari- 
ties of C/(q, z\Q) to the integral of Eq. (3.8) is strongly 



suppressed by the presence of the Bose factor (exp{/3z} — 
1) _1 therein, since j3\(i\ 3> 1 in the strong-coupling limit. 
In this limit, the contribution of these singularities can 
then be neg lect ed, with the result that the value of the 
expression (3/7) is proportional to the bosonic density 
Pb-i by the very argument discussed in the previuos sub- 
section for point-like bosons. 

We conclude that, for a system of composite bosons 
in the "low-density" limit, the leading diagrams contain 
just one cycle, like the ones shown in FigJ9b. In analogy 
with the corresponding diagrams of Fig. |p for point-like 
bosons, we shall refer to these diagrams as the T-matrix 
approximation for composite bosons. 

That in the strong-coupling limit all diagrams of 
Fig. ^|b are of the same order (in the original fermionic 
parameter kpap - cf. Ref. 14) can be further verified as 
follows. Consider the elementary rung (drawn in Fig. [|c) 
in the diagrammatic representation of the T-matrix for 
composite bosons. We know that the interaction ver- 



tex u 2 scales as \fj,\ 



-3/2 



x a F in the strong-coupling 



limit, where the factor 3 / 2 stems from the behaviour 
of ii2 (cf. the comment after Eq. ( 2.33| )) and the fac- 



tor ap fr om the presence of two powers of D in the 
definition (p.6| ) of the rescaled interaction vertex u 2 (cf. 
also Eq. ( |3 . 5[ ) ) . By the argument made for point-like 
bosons at the end of the previous subsection, the cut- 
off on the four-momentum of the rung is provided by 
the interaction vertex since the two "bare" propagators 
for composite bosons run in the same direction. This 
four-momentum sum is then proportional to the cut- 
off of the interaction vertex, which we know from Ref. 



16 to be given by 
and recalling that |/z 



Putting all factors together 
cx ap in the strong-coupling 
limit, we obtain that the rung of Fig. ^c scales like 



-1/2 



I Ml 



-3/2 



x a, 



constant. This result implies that 



the order of the diagram in the parameter kpap is not 
altered by the insertion of an arbitrary number of these 
rungs. 

The diagrams of Fig. |^ have been drawn in terms 
of propagators and interaction vertices for composite 
bosons. When extrapolating these diagrams from the 
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bosonic to the fermionic (weak-coupling) limit through 
the crossover region, however, it will be useful to draw 
them also in the alternative representation in terms of 
the constituent fermions, via the correspondence rules of 
Fig. ^. In Appendix B we will examine in detail which 
fermionic diagrams are associated with the T-matrix ap- 
proximation for composite bosons. In this way, we shall 
relate the value of the symmetry factors of the bosonic 
diagrammatic structure with the number of independent 
diagrams in the corresponding fermionic diagrammatic 
structure. We shall also show that in the weak-coupling 
limit, where the composite-boson propagator becomes a 
constant to leading-order in kpap (cf. Eq. ( 2.10D ), our 
generalized T-matrix approximation reduces correctly to 
the Galitskii approximation for the self-energy of a "low- 
density" Fermi gas. 



C. Coupled equations denning the generalized 
T-matrix approximation for the BE-BCS crossover 

Once the self-energy diagrams for composite bosons 
in the "low-density" limit have been selected according 
to the above prescriptions, there remains to determine 
the analytic expression of these diagrams. To this end, 
it is convenient to stick with the bosonic representation 
and write down the expression for the bosonic propaga- 
tor < b*(q)b(q) >5 ctf with the self-energy insertions of 
Fig. [|b, making use of the standard bosonic diagram- 
matic rules.E3'EJ Recalling the general relation (A7) of 
Appendix A for the "full" bosonic propagator (as well as 
its counterpart ( 2.26 ) for the "bare" bosonic propagator), 
we obtain: 



S 2 \. 

Si2u 2 (q', q, q, cf) - — ^ 2u 2 (q' , q, q%, q + q' - qi) 

91 

xY {q 1 )T {q + q' - qi)2u 2 (q + q' - qi, qi, q, q') 

+ (g V \2 ^2 2U2 (l'' 1' ?2, 9 + l' ~ 92)r o (g 2 )r o (g + q' - q 2 ) 

x2u 2 (q + q' - q 2 ,q 2 , 9i, q + q' - q{]T {q{) T (q + q' - q x ) 



x2u 2 (q + q' -qi,qi,q, q') 



(3.10) 



In this expression, Sl = 2 1 ~ L (L = 1, 2, • • •) is the sym- 
metry factor associated with the bosonic self-energy dia- 
gram of Fig. |^b containing L bosonic interaction vertices 
[cf . Appendix B] , u 2 is proportional to the effective two- 
boson interaction of Eq. ( |2.28 ) 

1 

pV 



u 2 (qi 



s(-k)e(k + q 2 )e(-k + qi - qi)e(k + q A ) ' 



(3.11) 



and r o is the "bare" particle-particle ladder of Fig. |^a, 
whose complete expression 



r D (g) 



tanh(/3f(k)/2) 



dk 

tanh(/3£(k- q)/2) m 

k2 



2(£(k)+£(k-q)-»n i ,) 



(3.12) 



holds irrespective of the value of clf and reduces in par- 
ticular to the form (3.4) in the strong-coupling limit 
j3fj, — > — oo. 

Note that the factor 2, which entered the definition 



( 2.28 ) for u 2 but has been removed from the defin ition 
( 3.11 ) for u 2l multiplies now each u 2 in Eq. ( 3.10 ) ex- 
plicitly. These fa ctors of 2 combine with the symmetry 
factors Sl in Eq. ( |3.10 ), resulting in an overall factor of 2 
multiplying the whole expression within brackets therein. 
This expression corresponds then to the following integral 
equation 



*b (91,92,93,94) = ^2(91,92,93,94) 
1 \ - 



PV ^ 

95 

xi B (qi 



02(91, 92, 95, 9i + 92 - 95^0(95) r o (9i + 92 - 95) 



92 - 95,95,93,94) 



(3.13) 



with q± + q 2 = (73 + 54, pr ovide d we set q\ = q± = q' and 
92 = 93 = 9- Equation ( |3.13| ) identifies the generalized 
T-matrix for composite bosons™ in . . an alogy with the ordi- 
nary T-matrix for true bosons BEJEII It is thus clear from 
Eq. ( 3.10p that the quantity 



^ ) (ff) = -^Ur„(,/)/ i> .«/.9-'/.'/) i^-N) 



represents the T-matrix approximation to the self-energy 
for composite bosons, in terms of which Eq. ( 3.10| ) ac- 
quires the characteristic form of a Dyson's equation: 



r(g)=r (g) + r o (q)X ( *\q)r o (q) + 
To(q) _ 1 



1 - Z%Hq)T (q) 



To (9) 



- 4%) 



(3.15) 



With the above expression for the particle-particle lad- 
der (which, by construction, correctly describes a system 
of "low-density" composite bosons in the strong-coupling 
limit), we can obtain the fermionic self-energy of inter- 
est in analogy to Eq. (2.35), by joining the incoming and 
outgoing arrows of the particle-particle ladder T with a 
single-particle fermionic Green's function in the only pos- 
sible way. We then write: 



(3.16) 



with r given by Eq. ( 3.15 ) and where Q° is the 
"bare" single-particle fermionic Green's function. The 
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self-energy (3.16) has in turn to be inserted into the 
fermionic Dyson's equation, to yield the full single- 
particle fermionic Green's function Q. Eventual ex- 
trapolation from the strong- to the weak-coupling limit 
through the crossover region requires us to eliminate the 
chemical potential in favor of the particle density p, by 
evaluating 



'5(k) 



(3.17) 



yond first order in the interaction potential for composite 
bosons, give contributions of the same order of magni- 
tude as the first-order (Hartee-Fock) diagram. In addi- 
tion, it will turn out that the series of diagrams depicted 
in Fig. does not converge, making any truncation of 
the series not appropriate. For this reason, it is essential 
to solve the complete integral equation associated with 
this series of diagrams to get a correct description of the 
strong-coupling limit. 



where rj is a positive infinitesimal. 

Besides the explicit Q° in Eq. ( 3.16| ), also all single- 
parti cle fe rmionic Green's functions entering the expres- 
sion ( 3.15 ) for r are meant to be "bare" ones, in anal- 
ogy to the original approach by Galitskii.EI We ex- 
pect, in fact, that the inclusion of self-consistency in 
the exp licit single-particle fermionic Green's function 
of Eq. ( 3.16| ) should not be essential to represent cor- 
rectly the fermionic self-energy, either in the strong- 
coupling limit (where, on physical grounds, it is rather 
the bosonic propagator that needs to be represented 
correctly) or in the weak-coupling limit (where self- 
consistency drops out anyway for a "low-density" Fermi 
system). [In addition, in the crossover region (where our 
"low-density" approximation merely interpolates between 
the two controlled limits) inclusion of self-consistency 
should not be of primary concern.] By the same token, in- 
clusion of self-consistency in the single-parti cle fer mionic 
Green's functions entering the expression ( 3.1 5| ) for T 
would yield contributions of higher order in the small 
parameter kpciF with respect to the ones retained (both 
in the weak- and strong-coupling limits). 

The question of self-consistency in the single-particle 
fermionic Green's function is also connected with the gen- 
eral requirement for an approximation to the fermionic 
self-enecgy to be "conserving" in the Baym-Kadanoff 
sense£3c3 According to this requirement, once an ap- 
proximation has been selected by physical considerations 
at the level of the fermionic two-particle Green's function 
(in the particle-particle channel), one should in princi- 
ple construct a functional $ such that the fermionic self- 
energy itself is, in turn, obtained by taking the functional 
derivative of $ with respect to the self-consistent Q. It is 
obvious that this procedure yields several additional self- 
energy terms, over and above the one obtained by joining 
the incoming and outgoing arrows of the particle-particle 
ladder with a single-particle fermionic Green's function, 
as we did in Eq. ( |3.16 ). It is, however, also clear that, by 
construction , these additional terms would be of higher- 
order in the "low-density" parameter, and can thus be 
consistently neglected. 

A complete numerical evaluation of Eqs. ( |3.11 )-( 3.17 ) 
exceeds the purposes of the present paper. In the next 
Section we shall calculate the scattering length for com- 
posite bosons in the strong-coupling lim it, as a degener- 
ate case of the T-matrix given by Eq. ( |3 .13 ). This cal- 
culation will ensure us that the diagrams of Fig. p|b, be- 



IV. NUMERICAL RESULTS FOR THE 
SCATTERING LENGTH OF COMPOSITE 
BOSONS 

The scattering length a characterizes the low-enersy. 
collisions for the scattering from an ordinary potential.ta 
For the mutual scattering of two particles (each of mass 
M), a can be expressed by the relation t(0) = 4na/M in 
terms of the ordinary T-matrix t(0) in the limit of van- 
ishing wave vector. In particular, within the Born ap- 
proximation t(0) is replaced by the Fourier trasform u(0) 
of the inter-particle potential for vanishing wave vector. 

In a similar way, we shall define the scattering length 
a b for composite bosons (each of mass 2m) in the strong- 
coupling limit and for vanishing density, by setting 



*b(0) 



4-7T 

2m 



(>B 



(4.1) 



(8n/(m 2 a F )) 2 i B (0) 



and t B (0) = 
Eqs. ( ^29| ), ( gggfl , and (jTlj|)]. This 



where t B (0) = 

tg (0,0, 0,0) [Cf. 

quantity is expected to be important for the calculation of 
the bosonic self-energy ( B.14| ) , insofar as the generalized 
T-matrix for composite bosons therein depends weakly 
on its arguments. 

To lowest order in the effective interaction for compos- 
ite bosons, we can replace t B (0) by 1(2(0) and write 



4-7T 

" 2 (0) = — aT n 
2m 



(4.2) 



withi n th e Born approximation. Comparison with 
Eq. ( 2.3C ) yields then the value a^ orn = 2ap, as an- 
ticipated by Eq. (2.31). Note also that the replacement 



of is by it 2 (or, equivalently, of i B by u 2 ) corresponds to 
the bosonic Hartree-Fock approximation, whereby only 
the diagram at the left in Fig. ^|b is considered instead of 
the whole sequence. 

To calculate i B (0), it is convenient to determine first 
the auxiliary quantity is (9, — q, 0, 0) by solving the fol- 
lowing closed-form equation 



t B (q, -0,0,0) = u 2 (q, -9,0,0) 



(4.3) 



E u 2 {q,-q,q\-q')T { q ')T {-q')i B { q ' ,-q' ,0,0) , 



which is obtained from Eq. (3.13) by setting q\ = —q 2 = q 
and (73 = 54 = 0. This integral equation can be solved by 
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standard numerical methods, e.g., by reverting it to the 
solution of a system of coupled linear equations. 

Before embarking into this numerical calculation, we 
can obtain a preliminary estimate of the value of is(0) 
with limited effort, by neglecting the four- vector depen- 
dence of Ib as well as the fr equ ency dependence of u 2 on 
the right-hand side of Eq. (4.3). We thus write approxi- 
mately: 



tfl(O) 
u 2 (Q) 

1 - 
1- 



t B (0) 

U2(0) 

-^3 ua(q, -q, 0, 0)- £ T (q)T (-q) 



dq 



(27T) 

dq 



3 u 2 (q,-q,0,0)-2^ 



^ v Am ^ v + Am 



■W2(q,-q,0,0)i 

2m 



(4.4) 



where the last result holds to the leading order in the 
density. Recalling further from Ref. 16 that 



«2(q, -q, 0,0) = w 2 (0) 



(IqM 



(4.5) 



in the strong-coupling limit, we obtain for the integral 
on the right-hand side of Eq. (4.4) the following value (in 
three dimensions) 



dq . 1 

—-3 U2 (q,-q,0,0) ^ 

V ' 2m 



u 3 (0) 



7TQ, p 



= 4 (4.6) 



where use has been made of the result (2.30). Compari- 
son with Eqs. (4.1) and ( fi~^ ) yields eventually 



as 

-.Born 



t B (P) 

ua(0) 



(4.7) 



which implies that the contribution of the series diagrams 
depicted in Fig. is of the same order of magnitude 
as the contribution of the first-order (Hartree-Fock) dia- 
gram therein. 

The above approximate calculation suggests us to con- 
sider in more detail the scattering problem for two true 
bosons (each of ma ss M ), mutually interacting via a po- 
tential of the form (|4.5|), namely, 



u(q) 



7 



87ra; 



M l + <7 2 /4 



(4.8) 



where q — \q\dF and 7 is a dimensionless c oupl ing con- 
stant (7 = 1 and M = 2m for the potential J4.5| )). Solv- 
ing for this simplified scattering problem will, in fact, be 
instructive to obtain the solution of the original scatter- 
ing problem (4.3) for composite bosons, since it will (i) 
suggest a nontrivial approximation to be carried over to 
the original equation (4.3), and (ii) assess whether the 
region of interest (7 ~ 1) belongs to the perturbative or 



nonperturbative regime of the scattering integral equa- 
tion. 

To this end, we recall the equation satisfied by the 
T-matrix for two-body scattering. In particular, it is 
sufficient to consider the following degenerate form 



i(q,0) = u{q) - 



dq' u(q-q')t(q',0) 



(27T) 



q f 

M 



(4.9) 



which resembles Eq. Qt.S| ) for composite bosons but lacks 
its dependence on Matsubara frequencies. As already 
noted, the scattering length is related to t(0) = i(0, 0) by 
the relation a = <(0)M/(47r). For a spherically symmet- 
ric potential, t(q, 0) = t(|q|,0) and the angular integral 
in Eq. (4.9) can be readily performed. The remaining 
integral over |q| can be suitably discretized over a mesh, 
until convergence is achieved for the desired value t(0) 
(200 mesh points have proved sufficient to get a 1% ac- 
curacy in the scattering length). 

The results for a (in units of oj?) vs. 7 are shown in 
Fig. where the asterisks corresponds to the numer- 
ical solution of Eq. fl4.9| ) and the full curve represents 
the Born approximation a Born /aF — 27. Note that, in 
the region of i nter est (7 as 1), a Boin /a as 3 and the so- 
lution of Eq. (O) cannot be obtained by perturbative 
methods (the perturbative region - where the Born se- 



ries associated with the integral equation (4.9) for a re- 
pulsive potential would converge - is, in fact, limited by 

a Born /a ~ 2 and corresponds to 7 ~ 0.5). For this rea- 
son, any truncation of the integral equation would not 
be justified. Note also that the difference between a Born 
and a increases drastically as 7 increases (we have veri- 
fied that a jap is proportional to log 7 at least over eight 
decades). 
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FIG. 10. Scattering length a (in uni ts of of) vs. the cou- 
pling strength 7 of the potential (|4.8|) , obtained by solving 
numerically the integr al equ ation ( |4.9| ) (asterisks) and by the 
approximate relation (4.11) (broken line). The full line rep- 
resents the Born approximation a/aF = 27. 



Reducing the integral equation to a matrix equation 
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(by discretizing the integral it contains) may, how ever , 
be less accurate when solving the original equation (4.3). 
In that case, it might be useful to make comparison with 
an alternative approximate result that avoids matrix in- 
version altogether. The solution of the two-body problem 
(4.9) suggests us to get this approximate result along the 
following lines. Quite generally, suppose that u(q) ha s 
a characteristic range q . It then follows from Eq. fl4.9| ) 
that t(q, 0) has the same range q and is proportional 
to u(q) for |q| 3> q (both properties can be verified by 
working out specific examples). Thus, if we set 



*(q,Q) 

t(0) 



A(q) 



u(0) 



(4.10) 



then A(q) varies slowly from A(0) = 1 to A(oo) < 1 (fo r 
a repulsive potential). Since the integral in Eq. (4.S) 
weighs more "small" than "large" wave vectors, we can 
replace A(q) ~ A(0) = 1 therein and obtain 



*(0) 



u 2 (0) 



'(0) + / 



2 (q) 



(4.11) 



(2ir) 3 q 2 /M 



The values of a/ of obtained with this approximation are 
reported in Fig. [l^ (broken line), and show remarkable 
agreement with the full numerical solution of the integral 
equation ( [4.9] ), from the perturbative region (7 <C 1) up 
to the region of interest (7 w 1). 

The same sort of approximation ca n be adopted when 
solving the original integral equation (4.3) for composite 
bosons. In this case, we set approximately 



i fl (g, - g ,0,Q) ua(?,-g,0,0) 



*b(0) 



«2(0) 



(4.12) 



and obtain from Eq. (4.3) the expression 

s 2 (0) 2 



*b(0) 



«2(0) + ua( ff , -q, 0, 0) 2 T (q) T (-q) 



(4.13) 



The calculatio n of U2(q, —q, 0, 0) according to the gen- 
eral expression (3.11) can be reduced to the numerical 
evaluation of an integral over |k|, with /i = — e a /2 in 
the fermionic single-particle Green's functions. The fre- 
quency part of the sum over q in Eq. (4.13) can further 
be reduced to an integral, since for vanishing density the 
range of temperature we are interested in vanishes. An 
accurate numeric al eva luation of the ensuing integral over 
|q| and f2 in Eq. (|4.13|) yields 



^Born 



*g(0) 
u 2 (0) 



1 



2.69 



(4.14) 



within an estimated 1% numerical accuracy. 

Full numerical calculation of Eq. (flj]) requires us to 
introduce a finite-size mesh for the variables (|q|, fi) as 



well as (|q'|, f2'), with the angular integral over q' affect- 
ing only the function U2(q,—q,q',~q')- Equation (4.3) 
is thus reduced to a set of coupled equations for the un- 
knowns tsflql, Q; |q|, — O; 0; 0), which we have solved by 
the standard Newton-Ralphson algorithm with a linear 
interpolation for the integral over |q'| and Q! . In this way 
we obtain 



-.Born 



t B (Q) 

«a(0) 



1 



2.65 



(4.15) 



within an estimated 5% n umeri cal ac curac y. The agree- 
ment between the results ( 4.14 ) and ( 4.15 ) implies that 
the numerical procedures we have adopted are accurate 
enough for the present calculation. 

To verify that this result could not be inferred from 
a perturbative expansion of the integral equation ( [4.3j ), 
we calculate eve ntu ally the second term on the right- 
hand side of Eq. fl4.3j ) by replacing ts(g', — q', 0, 0) therein 
with v,2(q', —q', 0, 0) and by setting q = eve rywh ere for 
convenience. In this way we obtain [cf. Eq. (|!l|: 



ifl(0) ~ ua(0) 



1 v ?2 2 (aW,0,0) 2 
PV Ga(0) 

q' 



Y (q')T {-q') + --- 
u 2 (0)(l - 1.69 + •••) 



(4.16) 



showing clearly that the geometric series would not con- 
verge in this case. 

To summarize, we have shown that, in the strong- 
coupling limit, the value as — 2a p obtained for 
the composite-boson scattering length within-, the self- 
consistent fermionic T-matrix approximation,!! is mod- 
ified to as — (3/4)ap by the correct inclusion of all low- 
density contributions for a system of composite bosons. 



V. CONCLUDING REMARKS 

In this paper, we have determined the correct diagram- 
matic approximation for a "dilute" system of composite 
bosons, which form as tightly bound pairs of fermions 
in the limit of strong attraction between the constituent 
fermions. We have emphasized that it is physically the 
comparison of the average interparticle distance to the 
characteristic length associated with the residual inter- 
action between the composite bosons to determine the 
"diluteness" condition in the strong-coupling limit of the 
original fermionic attraction. For this reason, the sim- 
ple argument, according to which to a "dilute" system 
of fermions in the weak-coupling limit there necessarily 
corresponds a "dilute" system of composite bosons in the 
strong-coupling limit, might be misleading, since two dif- 
ferent interactions actually control the "diluteness" con- 
dition in the two limits. It is thus essential to treat the 
residual interaction between the composite bosons with 
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the due care, in order to control the strong-coupling limit 
of the theory appropriately. In this context, it is worth 
mentioning that the importance of a proper treatment 
of the residual boson-boson interaction in the strong- 
coupling limit has been emphasized oil the pioneering 
paper by Nozieres and Schmitt-Rink,L3 but never duly 
taken into account in subsequent work. 

We have also shown that the selection of the diagram- 
matic contributions according to the "diluteness" pa- 
rameter proceeds along quite different lines in the weak- 
coupling limit (where the small parameter is kpap) and 
in the strong-coupling limit (where the small parameter 
is /5g 3 as). li-3 Accordingly, diagrammatic contributions of 

1/3 

the same order in p B og in the strong-coupling limit 
correspond, in general, to different powers of kpaF in 
the weak-coupling limit. Mathematically, this difference 
is due to the behaviour of the fermionic particle-particle 
ladder, which reduces to a constant in the weak-coupling 
limit but develops a singularity in the strong-coupling 
limit. 

Our selection of diagrammatic contributions has rested 
on a suitable regularization of the fermionic interaction, 
which has caused the ratio between the particle-particle 
and particle-hole contributions to be infinite. For a Hub- 
bard Hamiltonian on a lattice, where this regularization 
cannot be applied, we expect the difference between p- 
p and p-h contributions to be less extreme albeit still 
appreciable, so that our selectioB-pf diagrammatic con- 
tributions may still remain valid. E3 

Quite generally, we have remarked that, with our 
choice of the fermionic interaction, the most general 
structure of the diagrammatic theory is constructed with 
the "bare" particle-particle ladder plus an infinite set of 
(four, six, • • •, -point) vertices. This remains true for any 
value of the fermionic coupling and not just in the strong- 
coupling limit where the composite bosons form. We have 
also remarked that the "diluteness" parameter (kpap or 

1/3 

Pb °b) emerges naturally from the theory, both in the 
weak- and strong-coupling limits, without having to be 
imposed as an external condition. Accordingly, keeping 
tracks of the powers of this small parameter in the di- 
agrammatic theory can be relevant only in the weak- 
and strong-coupling limits. In the intermediate-coupling 
(crossover) region, on the other hand, a small parameter 
is lacking and consequently the diagrammatic approxi- 
mations cannot be controlled by any means .L3 

For these reasons, implementing the self-consistency of 
the fermionic Green's functions within the fermionic T- 
matrix approximation does not seem a priori to be an im- 
portant issue for the BCS-BE crossover. Self-consistency, 
in fact, drops out in the weak- and strong-coupling limits 
when the "diluteness" parameter is small, while in the in- 
termediate (crossover) region inclusion of self-consistency 
within the fermionic T-matrix approximation (as well as 
within any other approximation over and above it) can- 
not anyway be controlled by the lack of a small parameter 
(even though inclusion of self-consistency might produce 



in practice sizable numerical effects). 

We have emphasized in this paper that the (self- 
consistent) fermionic T-matrix approximation does not 
account properly for the boson-boson interaction in 
the strong-coupling limit, at least in three dimensions. 
This approximation, however, has been recently adopted 
to discuss pseudo-gap and related issues within, the 
negative-?/ Hubbard model in two dimensions oHd As- 
sessing to what extent the approach we have developed 
in this paper can be carried over to the two-dimensional 
case is not a priori evident and will require further in- 
vestigations. From physical intuition one would expect 
the bosonic regime to be reached even more effectively in 
two than in three dimensions, insofar as the two-fermion 
bound state is present in-two dimensions for any (attrac- 
tive) coupling strength.H3 Our dealing with the three- 
dimensional case first was required for manifesting at the 
outset the effects on the BCS-BE crossover due to the 
progressive formation of bound-fermion pairs, thus iso- 
lating them from other effects which are peculiar to the 
two-dimensional case. 

We have considered in this paper temperatures above 
T c only. Extending our approach below T c even for the 
three-dimensional case will also require further investiga- 
tions. In fact, divergencies occur in the particle-particle 
channel when approaching the critical temperature,EZI 
over and above the divergencies due to the formation of a 
bound state in the associated two-body problem, which 
we have taken care of in this paper. 

It is interesting to point out that, when written in 
terms of the constituent fermions (cf. Fig. [l3] below), 
the diagrammatic structure introduced by our approach 
bears some analogy with the parquet approximation uti- 
lized for the positive-?/ Hubbard model.El In addition, 
the level of numer ical effort required to solve our inte- 
gral equation Eq. ( 3.13| ) for the generalized T-matrix for 
composite bosons appears to be comparable with the ef- 
fort required to .solve the integral equations defining the 
parquet schemeO 

It is also interesting to point out the strong anal- 
ogy between the present treatment of the BCS to BE 
crossover in a condensed matter system and the so- 
called OAI mapoing introduced some time ago in nu- 
clear physicsJlj'Cj where a systematic mapping between 
the diagrammatic theories for (composite) bosons and 
(constituent) fermions is also provided, albeit in a quite 
different physical context and with the use of approxima- 
tions more specific to the nuclear problem. 

A final comment is in order about the emphasis we 
have given to the correct description of the boson-boson 
interaction in the strong-coupling limit. We know, in 
fact, that in the extreme strong-coupling limit the boson- 
boson interaction can be neglected and the composite 
bosons become effectively free. Since this extreme limit 
of free composite bosons is also retrieved by the (self- 
consistent) fermionic T-matrix approximation, it would 
appear that the improvement brought about by our the- 
ory (which takes specifically into account the interaction 
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between the composite bosons) affects oniy the way this 
extreme limit is approached. That the important physics 
may reside in the correct description of the approach to 
the asymptotic regimes and not merely of the asymp- 
totic regimes themselves, can be confirmed from previ. 
ous experience with the so-called Hubbard-III solution;c3 
in that case, neither the weakly-interacting (Fermi liq- 
uid) nor the strongly-interacting (Kondo) physics is re- 
produced, although the extreme (non-interacting and 
atomic) limits are both retrieved. 
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APPENDIX A: COMPOSITE-BOSON 
PROPAGATOR AND TWO-FERMION GREEN'S 
FUNCTION 

It was shown in Ref. 16 (cf. Appendix A therein) that 
the average va lue o f b(q — 0) (taken with the effective 
bosonic action ( 2 . 1 9|) ) is proportional to the average value 
of the operator (2.16) (taken with the original fermionic 
action ( 2.1 5| ) ) . Purpose of this Appendix is to show that 
a similar relation holds between the "full" bosonic prop- 
agator < b*(q)b(q) >s eH and the average value of the 
product of two operators (2.16). 

To this end, i t is c onvenient to add to the original 
fermionic action (2.15) the following (source) term 



rP 

5S = / drV[j(q,r)%T) + J*(q,r)B(q,r)] 

(Al) 

and to define the fermionic generating functional [cf. 
Eq. (p4j)] 



Z\J,J*} = 



J VcDc exp{-S- SS} 
J VcDc exp{-5} 



(A2) 



I ntrod ucing the Hubbard-Stratonovich transformation 
( 2.17 ), shifting the bosonic integration variables by let- 
ting b(q, r) — > 6(q, r) + J(q, r), and integrating out the 
Grassmann variables, we obtain the exact identity: 



SZ\J, J* 



5J(q,T')SJ*(q,T) 



<6*(q,r')fe(q,r) >< 



1 



(A3) 



Upon differentiating the original expre ssio n (A2) in a 
similar way and comparing with Eq. (A3), we obtain 
eventually the desired relation: 

< &*(q,r')6(q,r) > Srff = V 2 < B(q, r')B(q, r) > s 

-V6(t-t'), (A4) 

where the last term on the right-hand side can be 
dropped since it vanishes with the regularization (2.4) 
in the limit k a — > oo. 



The result (A4) can be conveniently interpreted in 
diagrammatic terms as follows. Takin g th e Matsubara 
Fourier transform of both sides of Eq. ( A4 ) yields 



< b*(q)b(q) > Sea = V 2 ™ ( k + 0/2) « (k' + q/2) 



k.k' 



x < c T (fc + q)c l (-k)c l (-k')c ] ik' + q) > S 

= V 2 Y, w ( k + q/ 2 ) w (k' + q/2) G 2 {k, k'; q) (A5) 



k,k' 



which defines the fermionic two-particle Green's func- 
tion Qi in the particle-particle channel. This function 
is depicted in Fig. |ll], in terms of the single-particle 
fermionic Green's functions Q and of the generalized 
particle-particle ladder V. Note that the external single- 
particle lines of Fig. [ll] (as well as the single-particle lines 
contained in T) include, in principle, self-energy inser- 
tions. 
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FIG. 11. Fermionic two-particle Green's function Q2, ex- 
pressed in terms of the "full" single-particle Green's function 
Q (represented here by a line) and of the generalized par- 
ticle-particle ladder F (represented by a box). The factor 
{f3V)~ originates from the diagrammatic rules. Spin labels 
are indicated by up and down arrows. 

Owing to our choice of a "contact" potential, the gen- 
eralized particle-particle ladder of Fig. [ll] depends only 
on the sum q of the incoming (outgoing) four-momenta. 
For this reason, the sums over k and k' in Eq. ( A5) are de- 
coupled and affect only the external single-particle lines. 
Let 

^p(«) = \ E w ( k + q/ 2 ) £( fc + «)£(- fc ) ( A6 ) 
p k 

be t he a nalogue of the quantity I PP (q) introduced in 
Eq. (2.5), where now the full single-particle fermionic 
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Green's function Q replaces its bare counterpart Q° . We 
know already that I pp = — 1 + 0(k^ 1 ), whereby only the 
behavior of Q° for large wave vectors contributes to the 
finite term. The self-energy insertions contained in Q can- 
not modify this result, since they vanish for large wave 
vectors (provided a well-defined finite, expression is asso- 
ciated with the self-energy itself in the limit k — > oo). 
This remark implies that also T pp = — 1 + 0(k~ 1 ). 

Combining Eq. ( |A5| ) with the representation of Fig. [Il], 
we then obtain the following relation between the "full" 
bosonic propagator and the generalized particle-particle 
ladder: 



< b*(q)b(q) >< 



(A7) 



of which Eq. (2.26) of the text is the simplest approxima- 
tion. It is clear from our derivation that the result (A7) 
holds irrespective of the value of the fermionic scattering 
length ap. 



The equivalence ( A7) between the generalized particle- 
particle ladder and the "full" bosonic propagator is basic 
to the purposes of the present paper, since it enables us 
to interpret in a rigorous fashion the general diagram- 
matic structure of T(q) in terms of the building blocks of 
the diagrammatic structure of < b*(q)b(q) >s etf , namely, 
the "bare" bosonic propagator < b*(q)b(q) >s 2 tt given 

-point interaction 
■, in the order. 



by Eq. (2.26) and the four , s ix, ■ ■ 
vertices given by Eqs. fl2 28[ ), fl2.33p , 



APPENDIX B: FERMIONIC DIAGRAMS 
CORRESPONDING TO THE T-MATRIX 
APPROXIMATION FOR COMPOSITE BOSONS 

In this Appendix we discuss the mapping between the 
diagrammatic structures of the composite bosons and the 
constituent fermions, for the case when only the four- 
point interaction vertex for composite bosons is retained 
and all higher-order vertices are neglected. In particular, 
we shall consider the class of diagrams of interest depicted 
in Fig. ^b, corresponding to the T-matrix approximation 
for composite bosons. To make our argument more com- 
plete, we shall also consider one additional bosonic di- 
agram not belonging to this class, owing to its peculiar 
topological structure. By working out these two exam- 
ples in detail, we shall establish a definite correspondence 
between the value of the symmetry factors for the bosonic 
diagrammatic structure and the number of independent 
diagrams in the associated fermionic structure, as an- 
ticipated in the Introduction. Besides providing a com- 
pelling check on the general structure of our mapping, 
determining the explicit form of the fermionic diagrams 
proves useful to verify that our generalized T-matrix ap- 
proximation correctly reduces to the standard fermionic 
(Galitskii) approximation in the weak-coupling limit. 

We begin by determining the symmetry factors as- 
sociated with the diagrams of Fig. for true bosons, 



where the two-body interaction is assumed to be suitably 
symmetrized. The symmetry factor iS^ associated with 
the diagram containing L interaction vertices is given by 
the ratio of the number of ways to construct this dia- 
gram from the corresponding prediagram and the factor 
A L L\ o*iginating from the perturbative expansion at or- 
der L.tSEy By definition, in the prediagram only the ex- 
ternal (incoming and outgoing) arrows of the self-energy 
and the interaction vertices appear, but not the propaga- 
tors that join the vertices among themselves or with the 
external arrows. Note that, for a complex bosonic field, 
the interaction vertices bear directional arrows and, con- 
sequently, the number of ways to join them with the prop- 
agators differs from the value of the standard <j) 4 theory 
with a real field. One can readily show that the num- 
ber of ways to construct the diagram of Fig. ||b with L 
interaction vertices from the corresponding prediagram 
is given by 2 L+1 Ll, yielding for the symmetry factor the 
value 



>£+i 



1 



4 L Ll 



2 L-l 



(Bl) 



Since for this set of diagrams all arrangements of vertices 
and propagators produce different contractions (accord- 
ing to Wick's theorem), no other factotis required in this 
case beside the symmetry factor (|B1| ).eZ1 

It is clear that with the diagrams of Fig. |^b for compos- 
ite bosons is also associated the same symmetry factor 
(Bl), as they share the same topological structure of the 
diagrams of Fig. Since the symmetry facte* of any 
fermionic diagram is instead bound to be unityj^j we ex- 
pect iS^ 1 identical fermionic diagrams to be associated 
with a given diagram for composite bosons having sym- 
metry factor Sl- 

To verify this statement, we recall the correspondence 
(shown in Fig. ||) of the propagator and vertex of the 
bosonic theory with the building blocks of the fermionic 
diagrammatic structure (which takes also into account 
the spin labels). In subsection IIB we have already dis- 
cussed the case with L = 1 and Si = 1, with which two 
distinct fermionic diagrams have consistently been asso- 
ciated (cf. Fig. |). The case with L = 2 and 6>2 = 1 /2 is 
shown graphically in Fig. [l^. In this case, it can be read- 
ily verified that the four possible fermionic diagrams (ob- 
tained by associating two distinct fermionic structures to 
each bosonic vertex via the rules of Fig. [|b) are equal in 
pairs. The resulting factor of 2 multiplying each indepen- 
dent fermionic diagram is then canceled by the symmetry 
factor 1S2 = 1/2, so that each independent fermionic dia- 
gram is correctly multiplied by unity, as anticipated. 

By a similar token, when L = 3 one can verify that, 
out of the eight possible fermionic diagrams, only two 
are independent. The ensuing factor of 4 in front of each 
of these two independent diagrams is then canceled by 
the symmetry factor S3 = 1/4. More generally, if we 
call "up" and "down", respectively, the two fermionic 
structures associated in Fig. |4|b with the bosonic four- 
point vertex, it can be verified that fermionic diagrams, 
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obtained from each other by replacing a pair of "up" 
by a pair of "down" fermionic structures, are equal. 
This property is sufficient to guarantee that each dis- 
tinct fermionic diagram is correctly multiplied by unity. 
Consider, for instance, two nontrivial cases with odd 
(L = 5) and even [L = 6) number of bosonic vertices. 
When L = 5, 2 5 = 32 fermionic diagrams are gen- 
erated from the bosonic diagram. In this case, from 
the single fermionic diagram with all "up" structures 
we can generate 5!/(3!2!)=10 diagrams with three "up" 
and two "down" structures, and 5!/(l!4!)=5 diagrams 
with one "up" and four "down" structures, for a total 
of 1+10+5=16 identical fermionic diagrams. The same 
number of identical diagrams results from the other in- 
dependent diagram having all "down" structures (which, 
evidently, cannot be obtained from the previous diagram 
with all "up" structures by replacing an even number 
of "up" structures by "down" structures). The resulting 
factor of 16 associated with each of these two independent 
fermionic diagrams is then canceled by the symmetry fac- 
tor £5 = 1/16, and each independent fermionic diagram 
is correctly multiplied by unity. 



independent diagrams (out of the 64 possible fermionic 
diagrams) have either all "up" structures or five "up" 
and one "down" structures. To the first case we as- 
sociate l+6!/(4!2!)+6!/(2!4!)+l=32 identical diagrams, 
and to the second case 6!/(5!l!)+6!/(3!3!)+6!/(l!5!)=32 
identical diagrams. The resulting factor of 32 multiply- 
ing each independent fermionic diagram is then canceled 
by the symmetry factor Sq = 1/32, and each independent 
fermionic diagram is again correctly multiplied by unity. 

Quite generally, one can show that the fermionic dia- 
grams associated with the T-matrix diagrams of Fig. 
for composite bosons can be organized into three classes, 
containing zero, one, and two fermionic loops, respec- 
tively. A representative diag ram of each of these three 
classes is depicted in Fig. [l3pi, Fig. |l3|b, and Fig. |lg|c, in 
the order. Note that fermionic diagrams containing zero 
or two loops are associated with bosonic diagrams with 
an even number of interaction vertices, while fermionic 
diagrams containing one loop are associated with bosonic 
diagrams with an odd number of interaction vertices (cf. 
the comment made in subsection IIB while discussing 
Fig. |). 
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FIG. 12. The two distinct fermionic diagrams correspond- 
ing to the bosonic diagram with two interaction vertices 
(L = 2 and <S 2 = 1/2) depicted in Fig. §d. Boxes and full 
lines stand for the "bare" particle-particle ladder T and for 
the "bare" single-particle fermionic Green's function Q°, re- 
spectively. 

When L = 6, 2 6 = 64 fermionic diagrams are gener- 
ated from the bosonic diagram. In this case, the two 



(C) 



FIG. 13. Representative diagrams of the three classes of 
fermionic diagrams associated with the bosonic T-matrix, 
containing (a) zero, (b) o ne, and (c) two fermionic loops. 
Conventions are as in Fig. [L2l 

The above analysis can be made more systematic, by 
allowing also for the possibility of "twisting" pairs of 
fermionic lines while connecting the fermionic particle- 
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particle ladder of Fig. [|a with the four-point vertex of 
Fig. |]b. In other words, one may consider the situation 
where two pairs of fermionic lines, which would not be 
allowed to connect with each other "directly" as their 
spins would not match, are instead allowed to connect 
by "twisting" one pair of lines with respect to the other 
one. A careful analysis of the unraveling of the fermionic 
diagrams containing "twistings" , however, shows that no 
new diagram is introduced in this way, over and above the 
diagrams depicted in Fig. [l3] which have been generated 
by not allowing "twisted" connections. 

We are now in a position to verify that the fermionic 
diagrams, associated with the generalized T- matrix ap- 
proximation for composite bosons, correctly reduce to the 
standard Galitskii approximation in the weak-coupling 
limit, in the sense that they yield contributions of pro- 
gressively higher order in the small parameter kpaF with 
respect to the "bare" particle-particle ladder. For in- 
stance, while the particle-particle ladder is proportional 
to of in the weak-coupling limit (cf. Eq. ( 2.10p ), correc- 
tions to this ladder due to the simple diagrams of Fig. | 
are readily seen to be smaller by a factor (kpap) 2 . Sim- 
ilarly, corrections to the particle-particle ladder due to 
the diagrams of Fig. [12] are smaller by a factor {kpap) A 
with respect to the "bare" ladder. In general, to esti- 
mate in the weak-coupling limit the order in kpap of a 
given diagram contributing to the generalized particle- 
particle ladder T, we can rely on dimensional considera- 
tions , by counting the powers of ap in terms of the num- 
ber of "bare" ladders T a and the powers of kp in terms 
of the dimensionality of the four-vector sums over prod- 
ucts of internal single-particle fermionic Green's func- 
tions (which can safely be done because all sums are 
convergent). By this procedure, we obtain that, in the 
weak-coupling limit, the diagram of order L within the 
bosonic T-matrix approximation is smaller by a factor 
(kpap) 2L with respect to the "bare" ladder r o . 

This proves that diagrams of the same order in p 1 J 3 as 
in the bosonic (strong-coupling) limit end up having dif- 
ferent orders in kpap in the fermionic (weak-coupling) 
limit; accordingly, they would have been dismissed as 
being irrelevant, if the selection of diagrams would have 
been made directly for the weak-coupling limit. 

Finally, let us consider the bosonic diagram of Fig. |lj 
which does not belong to the set of T-matrix diagrams. 
For this diagram, a naive counting of the number of ways 
it can be constructed from the corresponding prediagram 
would give unity for the symmetry factor. However, a 
more careful analysis shows that a permutation of the 
two internal vertices is actually equivalent to a renaming 
the propagators joining these vertices, and has thus not 
to be regarded as an independent permutation.cll This 
implies that the bosonic symmetry factor for the diagram 
of Fig. |l4| is 1/2 and not 1. As a matter of fact, when 
generating the fermionic diagrams associated with this 
bosonic diagram according to the correspondence rules 
of Fig. ||, one finds that the 2 4 = 16 diagrams (produced 



by the presence of four bosonic interaction vertices) are 
equal in pairs, yielding just 8 distinct fermionic diagrams 
in agreement with the value 1/2 of the bosonic symmetry 
factor. As before, this result is not affected by introduc- 
ing "twistings" . 

Although we have examined in detail only two kinds of 
bosonic diagrams, we expect the correspondence between 
the bosonic symmetry factor and the number of indepen- 
dent fermionic diagrams to remain valid for all possible 
diagrams. 




FIG. 14. Bosonic diagram with four interaction vertices, 
which does not belong to the set of T-matrix diagrams and 
possesses a special internal symmetry. 
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